第 40 卷 第 2 其 应 用 力 学 学 报 .Yol40No.2 
2023 年 4 月 Chinese Journal of Applied Mechanics “一 [IC ND | FA APr. 2023 


文章 编号 :1000- 4939(2023 )02- 0310- 10 


非 比 例 阻尼 体系 等 效 Rayleigh 阻尼 
模型 时 域 计 算 方法 
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摘 要 :Rayleigh 阻尼 模型 具有 数学 简易 性 的 优点 ,应 用 广泛 ,其 阻尼 短 阵 构造 依赖 于 结构 模 态 阻尼 
比 。 结 构 阻 尼 ( 复 阻尼 ) 模 型 的 阻尼 短 阵 直接 由 材料 损耗 因子 和 刚度 矩阵 决定 ,在 非 比例 阻尼 体系 
中 具有 阻尼 矩阵 便于 构造 的 优点 ,但 存在 时 域 计 算 结 果 发 散 、 初 值 条 件 不 易 确 定 和 频 响 函数 非 因 果 
等 缺陷 。 本 研究 结合 两 种 阻尼 模型 的 优点 ,分 别 依据 阻尼 衰减 和 阻尼 耗 能 ,提出 了 与 结构 阻尼 模型 
,等 效 的 Rayleigh 阻尼 模型 。 算 例 分 析 结 果 表 明 : 等 效 Rayleigh 阻尼 模型 克服 了 结构 阻尼 模型 的 缺 
“ 陷 , 同 时 保留 了 非 比 例 阻尼 体系 中 阻尼 矩阵 易 构 造 的 优点 ;与 基于 阻尼 衰减 等 效 的 Rayleigh 阻尼 模 
型 相 比 ,基于 阻尼 耗 能 等 效 的 Rayleigh 阻尼 模型 计算 结果 近似 相等 ,但 避免 了 复 模 态 分 析 , 且 计算 
过 程 直 观 简单 。 
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Time-domain calculation method of equivalent Rayleigh damping 
model for non-proportional damping system 
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Abstract:Rayleigh damping model has the advantage of mathematical simplicity, but the construction of 
damping matrix depends on the the structural modal damping ratios. The damping matrix of the structural 
damping model (complex damping model) can be directly constructed by material loss factors and struc- 
tural stiffness matrix ,which has the advantage that the damping matrix is easily constructed for non-propor- 
tional damping system. However ,the structural damping model has some shortcomings , such as the diver- 
gence of the time-domain calculation resujlts ,the difficulty of determining the initial value conditions and 


the non-causality of the frequency response function. Combined with the advantages of the two damping 
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models ,the Rayleigh damping model is equivalent to structural damping ,which is constructed in this paper 


according to the equivalent damping attenuation and equivalent damping energy dissipation ,respectively. 


The numberical examples show that the equivalent Rayleigh damping model can overcome the shortcomings 


of the structural damping model and retain the advantage that the damping matrix can be easily constructed 


for the non-proportional damping system. Compared with the Rayleigh damping model based on the equiva- 


lent damping attenuation, the calculated results of the Rayleigh damping model based on the equivalent 


damping energy dissipation are approximately equal. Moreover,complex modal analysis is avoided and the 


calculation process is more intuitive and simple. 


Key words :non-proportional damping ;structural damping ; Rayleigh damping; damping attenuation ; damp- 


ing energy dissipation 


随 着 结构 功能 性 的 增多 和 安全 性 的 提升 ,单一 
材料 组 成 的 比例 阻尼 结构 较 难 满足 人 们 对 工程 结构 
的 要 求 , 由 不 同 阻尼 特性 材料 组 成 的 非 比 例 阻尼 混 
合格 构 逐 渐 被 推广 使 用 "”。 因 此 , 非 比 例 阻尼 体系 
动 为 响应 的 计算 方法 研究 具有 重要 意义 。 灰 性 阻尼 
模 负 和 结构 阻尼 模型 中 (或 称 复 阻尼 模型 "1 、 滞 变 
随 志 模型 '9 、 率 无 关 阻尼 模型 "等 ) 是 最 常用 的 结 
梅内 阻 尼 模 型 , 黏 性 阻尼 模型 具有 计算 过 程 简单 
的 已 点 ,但 非 比例 阻尼 体系 的 阻尼 算 阵 不 再 满足 经 
上 典 息 交 条 件 ,需要 借助 于 状态 空间 法 ,采用 复 模 态 
三 加 法 求解 结构 动力 响应 !"] 。 同 时 ,基于 秋 性 阻尼 
模型 的 非 比 例 阻 尼 体 系 阻尼 和 矩阵 构造 十 分 困难 。 
与 籍 性 阻尼 模型 相 比 , 结 构 阻 尼 模型 具有 阻尼 矩阵 
易 移 造 的 优点 ,其 阻尼 矩阵 可 直接 依据 材料 损耗 因 
子 构造 得 到 ”1 ,避免 了 结构 模 态 阻尼 比 的 测定 问 
题 志 基于 结构 阻尼 模型 的 复 模 态 亚 加 法 可 直接 在 物 
理 盈 间 内 实现 ,但 结构 阻尼 模型 存在 时 域 发 
散 C(3 7 \ 频 响 函 数 非 因果 ”有 阻尼 自 振 频 率 随 
损耗 因子 增 大 而 增 大 5 等 缺陷 ,这 使 得 结构 阻尼 模 
型 仅 适 用 于 有 限 频 率 范围 内 结构 的 频 域 分 析 21 。 
Rayleigh 阻尼 模型 是 一 种 简单 .应 用 广泛 的 符 
性 阻尼 模型 ,通过 假定 阻尼 矩阵 与 质量 矩阵 .刚度 矩 
阵 成 比例 ,保留 了 其 数学 简易 性 的 优点 。 为 了 计算 
非 比 例 阻尼 混合 结构 的 动力 反应 ,常见 做 法 是 将 非 
比例 阻尼 体系 近似 等 效 为 比例 阻尼 体系 ,采用 Ray- 
leigh 阻尼 模型 求解 动力 响应 。 但 Rayleigh 阻尼 矩阵 
的 构造 不 仅 依赖 于 质量 矩阵 和 刚度 矩阵 ,还 取决 于 
结构 的 模 态 阻尼 比 。 对 于 大 型 复杂 结构 , 模 态 阻 尼 
比 的 测定 是 十 分 困难 的 ,通过 振动 试验 也 仅 能 得 到 
低 阶 振 型 的 模 态 阻尼 比 。 因 此 ,在 结构 阻尼 模型 的 


为 了 利用 结构 阻尼 模型 中 阻尼 和 矩阵 易 构 造 的 优 
点 ,文献 [22 ] 利用 振 型 向 量 构建 权重 函数 ,得 到 了 
等 效 于 结构 阻尼 模型 的 Rayleigh 阻尼 模型 ;文献 
[23 ] 以 振 型 参与 质量 系数 作为 权 系数 ,采用 最 小 二 
乘法 将 Rayleigh 阻尼 和 矩阵 与 结构 阻尼 和 矩阵 等 效 , 进 
而 得 到 等 效 Rayleigh 阻尼 模型 。 但 上 述 方法 的 等 效 
过 程 较为 复杂 ,不 适宜 在 工程 实践 中 推广 使 用 。 

针对 非 比例 阻尼 体系 的 动力 响应 求解 问题 ,本 
研究 在 结构 阻尼 模型 的 基础 上 ,依据 结构 瞬 态 反应 
和 稳 态 反应 各 自 的 特点 ,分 别提 出 基于 阻尼 衰减 等 
效 的 Rayleigh 阻尼 模型 和 基于 阻尼 耗 能 等 效 的 Ray- 
leigh 阻尼 模型 ,可 克服 结构 阻尼 模型 时 域 发 散 的 缺 
隐 , 同 时 保留 了 其 阻尼 矩阵 易 构 造 的 优点 ,是 一 种 效 
率 更 高 的 非 比 例 阻尼 体系 的 时 域 动力 响应 计算 方法 。 


1 结构 阻尼 模型 的 缺陷 分 析 


基于 结构 阻尼 模型 的 单 自由 度 体系 运动 方程 为 
mx(t) + kx(t) + imkx(t) =f+if (1) 
式 中 : m 为 结构 的 质量 ;% 为 结构 的 刚度 ; 7 为 损耗 
因子 ; /为 外 激励 ;了 为 /的 复 化 对 偶 项 ;i 为 虚数 
单位 ， 即 1 = Vv-1。 
x(1) 为 结构 位 移 响 应 的 复数 表达 形式 , 即 
X(t) = Xi(t) + iw,(t) (2) 
式 中 ,x(t) 为 结构 的 真实 位 移 响 应 。 
式 (1) 对 应 的 自由 振动 方程 为 
mx(t) + kx(t) +imkx(t) =0 (3) 
采用 复 平面 法 求解 式 (3) ,可 得 
xi(t) = [Acos(wi) +Asin(w.it) Je ”+ 


[Ascos(@w.t) + Asysin(w.t) je (4) 


基础 上 ,建立 合理 ,等 效 的 Rayleigh 阻尼 模型 能 有 效 
提高 非 比 例 阻 尼 体 系 动力 响应 计算 的 效率 。 


式 中 : 4 、4，、4; 和 4 为 待定 系数 ; a, 为 自由 衰减 
系数 ; w, 为 有 阻尼 自 振 频率 。 
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为 有 阻尼 上 自 振 频率 。 
a, = 《mw 
| (10) 
w, =w Vl- 
式 (9) 和 式 (10) 表 明 , 相 比 结构 阻尼 模型 ,条 性 阻尼 


由 式 (4) 和 式 (5) 可 知 ,结构 阻尼 模型 计算 的 自 
由 振动 通 解 中 包含 发 散 项 ,导致 结构 阻尼 模型 的 时 
域 计 算 结果 存在 发 散 现象 ,不 能 稳定 收敛 ; 式 (4) 中 
存在 4 个 待定 系数 ,而 动力 计算 的 初 值 仅 包含 初始 
位 移 x,(t) 和 初始 速度 x1(to) ,导致 结构 阻尼 模型 
ee < 件 不 易 确 定 。 式 (5 ) 表明 

于 结构 阻尼 模型 的 有 阻尼 自 振 频率 随 着 损耗 因子 
ee 
非 因果 的 缺陷 “-”!。 因 此 ,结构 阻尼 模型 无 法 直接 
用 二 计 算 结构 的 时 域 动力 响应 。 


2 二 基于 阻尼 衰减 等 效 的 Rayleigh 阻 
马 尼 模 型 


、“ 阻 尼 将 导致 结构 的 自由 振动 响应 逐渐 收敛 为 静 
此 状态 ,因此 ,本 研究 首先 从 自由 振动 衰减 出 发 , 构 
建 等 效 于 结构 阻尼 模型 的 黏 性 阻尼 模型 。 然 后 ,在 
蕊 各 由 度 体系 的 等 效 入 性 阻尼 模型 基础 上 ,依据 多 
所 放 度 体系 的 结构 阻尼 模型 ,进一步 构建 非 比 例 阻 
尼 体 系 的 等 效 Rayleigh 阻尼 模型 。 


2: 单 自由 度 体系 的 等 效 医 性 阻尼 模型 


由 式 (4) 可 知 ,舍弃 结构 阻尼 自由 振动 通 解 中 
nh 页 ,可 得 到 稳定 收敛 的 部 分 自由 振动 通 解 为 

wi(t) = [Aicos(w.t) + hsin(w.t) Je”” (6) 

基于 黏 性 阻尼 模型 的 单 自由 度 体系 自由 振动 方 
程 为 


mi + hx(t) + ei(t) =0 (7) 
式 中 ; m 为 结构 的 质量 ; 上 为 结构 的 刚度 ; 。 为 结构 
的 黏 性 阻尼 系数 , 且 
= (8) 
CU 


式 中 : 为 结构 的 阻尼 比 ; w 为 结构 的 无 阻尼 自 振 频 


率 , 且 w = ,/ 二 。 
m 


求解 式 (7) ,可 得 到 基于 黏 性 阻尼 模型 的 结构 
自由 振动 通 解 为 

x(1) = [Bicos(wi) +Psin(wi)je (9) 
式 中 : B, 和 B, 为 待定 系数 ; a, 为 自由 衰减 系数 ; w， 
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模型 的 自由 振动 通 解 中 仅 包含 指数 衰减 项 ,不 包含 
指数 增加 项 , 且 有 阻尼 自 振 频 率 随 阻尼 比 的 增 大 而 
减 小 。 为 保证 两 种 阻尼 模型 的 阻尼 衰减 相等 , 即 自 
由 衰减 系数 相等 ,由 式 (5) 和 式 (10) 可 得 


cmw = w [tml (11) 
进一步 得 到 
< = [tn 1 (12) 


由 式 (12) 可 得 等 效 犁 性 阻尼 运动 方程 为 
i =f (13) 


2.2 多 自由 度 体系 的 等 效 Rayleigh 阻尼 模型 
于 结构 阻尼 模型 的 多 自由 体系 时 域 运 动 方 


mxX(t) + kx(t) + 


程 为 

MxX(i) + Kx(it) +iK.x(t) =F+iF’ (14) 
式 中 :M 为 结构 的 质量 矩阵 ;天 为 结构 的 刚度 矩阵 ; 
K, 为 结构 的 阻尼 矩阵; 为 外 激励 向 量 ; 了 "为 的 
复 化 对 偶 项 。 

阻尼 矩阵 可 表示 为 


s 
K, = Dn k 


式 中 : 9。 为 第 i 种 材料 对 应 的 损耗 因子 ; K; 为 第 i 
种 材料 对 应 的 刚度 矩阵 。 
方程 式 (14) 对 应 的 特征 方程 为 


(15) 


MA +iK,+K=0 (16) 
求解 方程 式 (16) ,可 得 特征 根 为 
、 = 一 di 十 ib, (17) 
A 一 人 党 10) 


式 中 :特征 根 的 虚 部 5 为 混合 结构 第 j 阶 振 型 有 阻 
尼 上 自 振 频率 ;特征 根 的 实 部 aj 为 衰减 系数 。 

基于 黏 性 阻尼 模型 的 多 自由 度 体系 时 域 运动 方 
程 为 


MXx(1) + Kx(1) +Cx(t) =F (18) 
式 中 ,C 为 结构 的 阻尼 和 矩阵 。 
Rayleigh 阻尼 模型 具有 数学 表达 简单 等 优点 ， 


] 力 学 学 报 


第 2 期 


在 多 自由 度 结构 体系 的 动力 响应 计算 中 得 到 了 广泛 
使 用 。 阻 尼 和 矩阵 可 表示 为 


C=aM +BK (19) 
第 7 阶 振 型 对 应 的 衰减 系数 为 
oj = 全 + 下 (20) 


式 中 , wj 为 混合 结构 第 j 阶 振 型 无 阻尼 自 振 频 率 。 
令 Rayleigh 阻尼 模型 和 结构 阻尼 模型 的 第 j 阶 
振 型 自由 衰减 系数 相等 ,可 得 


CQ 1 2 
0 7 + oho; (21) 
选择 结构 某 两 阶 重要 振 型 的 自 振 频率 w， 和 

ww, , 式 (21) 可 进一步 转化 为 


1 2 
+ FBon = an 


人 
人 (22) 

求解 式 (22) ,可 得 到 Rayleigh 阻尼 和 矩阵 的 系数 
ac 条 为 


2 入 
2a, 0 二 2an,0, 


CU 2 CU 
m n 
(23) 
2an 加 2a， 
B 2 2 
Wn, 一 On 


CN 综 上 ,利用 自由 衰减 相等 ,可 构建 出 等 效 于 结构 
阻尼 模型 的 Rayleigh 阻尼 模型 。 


3 基于 阻尼 耗 能 等 效 的 Rayleigh 阻 
三 尼 模 型 


号 结构 的 动力 响应 是 由 有 瞬 态 响应 和 稳 态 响应 组 
成 ,本 研究 在 第 2 节 依据 瞬 态 响应 衰减 ( 自由 振动 响 
应 衰减 ) 构 建 了 与 结构 阻尼 模型 等 效 的 Rayleigh 阻 
尼 模型 。 本 节 则 根据 稳 态 响 应 的 特点 ,首先 从 阻尼 
耗 能 出 发 ,构建 等 效 于 结构 阻尼 模型 的 黏 性 阻尼 模 
型 。 然 后 ,在 单 自由 度 体系 的 等 效 黏 性 阻尼 模型 基 
础 上 ,依据 多 自由 度 体系 的 结构 阻尼 模型 ,进一步 构 
建 非 比 例 阻 尼 体 系 的 等 效 Rayleigh 阻尼 模型 。 


3.1 单 自 由 度 体系 的 等 效 黏 性 阻尼 模型 


基于 黏 性 阻尼 模型 的 单 自 由 度 体系 时 域 运 动 方 
程 为 


mx(t) + kx(t) + cx(t) =f (24) 
谐 波 作用 下 单 自由 度 体系 的 位 移 响 应 为 
x(t) = XsinOt (25) 
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式 中 :站 为 结构 位 移 响 应 的 振幅 ; 9 为 外 激励 谐 波 的 
谐 波 作用 下 ,基于 黏 性 阻尼 模型 的 单 自 由 度 体 
系 每 周期 耗 散 能 量 > 


2 


AW, = [WOENG (26) 
将 式 (25 ) 代 入 式 (26) ,可 得 
AW, = mcOX (27) 


基于 结构 阻尼 模型 的 单 自由 度 体系 时 域 运 动 方 
程 在 实数 域内 的 表达 式 为 


Wn + kxi(t) —- nkxs(t) =f (28) 
ms(t) + hrxs(t) + nkx(t) = 三 
谐 波 作用 下 单 自由 度 体系 的 位 移 响 应 为 
xi(t) = XsinQt 
| (29) 
xs(t) = — XecosOt 


谐 波 作用 下 ,基于 结构 阻尼 模型 的 单 自 由 度 体 
系 每 周期 耗 散 能 量 为 


AW, = -ia(Ddu() (30) 
将 式 (29) 代 入 式 (30) ,可 得 
AW, = mnkX (31) 


当 外 激励 频率 与 结构 的 自 振 频率 相等 (结构 共 
振 ) , 即 9 = w 时, 式 (27) 可 进一步 转化 为 


AW, = TCX (32) 
令 两 种 阻尼 模型 中 阻尼 耗 散 能 量 相等 ,可 得 
TcwX” = mnkX (33) 


进一步 得 到 
共有 三 2 (34 ) 
(CU 
由 式 (34) 可 得 等 效 黏 性 阻尼 运动 方程 为 
a + (1) aF 


3.2 多 自由 度 体系 的 等 效 Rayleigh 阻尼 模型 


基于 黏 性 阻尼 模型 的 多 自由 度 体 系 无 阻尼 振 型 
向 量 为 

了 =[P，P : 9Jv] (36) 
以 第 j 阶 振 型 自 振 频 率 为 振动 频率 的 谐振 作用 
下 ,多 自由 度 体系 的 位 移 响 应 与 振 型 形状 成 比例 , 即 
Xi) = YPsinot (37) 

式 中 ,P 为 位 移 啊 应 的 幅 值 系数 。 

结构 的 模 态 阻尼 耗 散 能 量 ， 
AW,; = Towi 户 PCp， (38 ) 
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结构 阻尼 模型 下 ,结构 的 模 态 阻尼 耗 散 能 量 为 
AW,, = TPP; Ky, (39) 
令 两 种 模型 的 模 态 阻尼 耗 散 能 量 相等 , 即 
1 
9 


站 二 
人 (40) 
J 
基于 Rayleigh 阻尼 模型 , 式 (40) 可 进一步 表示 为 
下 
KK 
ap Mp +ppKP = (41) 
J 


选择 两 个 结构 重要 的 自 振 频 率 w。 和 w,, 式 
(41) 可 进一步 转化 为 


hE 


K 
ap, Mo9,, 看 Bp, Koy, 一 一 全 
(42) 
oT K o 
ap, Moy, +BoKy, ee 
式 《42) 可 进一步 转化 为 
下 Dn 
QM +B 0 
(43) 
下 Zn 
十 = 一 
QM Bb n w, 
a 
oS 人 
CA = of 
| Pm Mp, (44) 
nL 一 op.Moy, 
k, = OK 
be (45) 
k, = pKy, 
v =o.K 
| m nm (46) 
he V, = 9 Ky, 
《 = 求解 式 (43 ) ,可 得 Rayleigh 阻尼 系数 为 
Kd 人 V 
—k, 于 —k 
二 Wi, Cn 
mk, mk 
(47) 
Dn Vn 
—m. 一 —m 
-OO cn 
B mh, 一 7 ke, 


因此 ,利用 模 态 阻尼 耗 散 能 量 相等 ,可 构建 出 与 
结构 阻尼 模型 等 效 的 Rayleigh 阻尼 模型 。 


4 算 例 分 析 


4.1 单 自 由 度 体系 算 例 
基于 阻尼 衰减 等 效 的 黏 性 阻尼 模型 ( viscous 


time-damain damping model based on free vibration re- 
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sponse, VFT) 与 结构 阻尼 模型 的 自由 衰减 系数 比 
值 为 


pr=1 (48) 

基于 阻尼 耗 能 等 效 的 黏 性 阻尼 模型 ( viscous 

time-damain based on the equivalent damping energy 

consumption ,VET) 与 结构 阻尼 模型 的 自由 衰减 系数 
比值 为 


7 49 

V2 V1i+n -2 人 

依据 式 (48) 和 式 (49 ) 可 计算 不 同等 效 阻尼 模 
型 的 自由 衰减 系数 比值 的 变化 规律 ,计算 结果 如 图 
1 所 示 。 随 着 损耗 因子 的 增加 ,基于 阻尼 衰减 等 效 
的 黏 性 阻尼 模型 与 结构 阻尼 模型 的 自由 衰减 系数 比 
值 恒 等 于 1 ,基于 阻尼 耗 能 等 效 的 黏 性 阻尼 模型 与 结 
构 阻 尼 模型 的 自由 衰减 系数 比值 接近 1。 因 此 ,两 种 
等 效 黏 性 阻尼 模型 的 自由 振动 响应 与 结构 阻尼 模型 
自由 振动 响应 中 剔除 发 散 项 的 计算 结果 近似 相等 。 


Pp2 = 


0.0 0.2 0.4 


0.6 
损耗 因子 
图 1 不 同等 效 阻尼 模型 的 衰减 系数 比值 变化 规律 


Fig.1 The variation laws of attenuation coefficient ratios 


0.8 1:0 


based on different equivalent damping models 
VFT 与 结构 阻尼 模型 的 动力 放大 系数 比值 为 
(1 - 8 ) 2 二 这 


人 位 8) +2( Vit -1)8 


(50) 
式 中 ,6 为 外 激励 谐 波 频 率 与 自 振 频率 的 比值 。 
VET 与 结构 阻尼 模型 的 动力 放大 系数 比值 为 
1 -6) +7 
ga 人 有 7 
依据 式 (50) 和 式 (51) 可 计算 不 同等 效 阻尼 模 
型 的 动力 放大 系数 比值 的 变化 规律 ,计算 结果 如 图 
2 所 示 。 随 着 损耗 因子 的 增加 ,基于 阻尼 衰减 等 效 
的 黏 性 阻尼 模型 与 结构 阻尼 模型 的 动力 放大 系数 比 
值 、 基 于 阻尼 耗 能 等 效 的 黏 性 阻尼 模型 与 结构 阻尼 
力学 学 报 


(51) 


模型 的 动力 放大 系数 比值 均 接 近 1。 因 此 ,两 种 等 
效 黏 性 阻尼 模型 的 稳 态 振动 响应 与 结构 阻尼 模型 稳 
态 振 动 响应 的 计算 结果 近似 相等 。 

地 震 作 用 下 结构 的 动力 响应 由 自由 振动 响应 和 
稳 态 振动 响应 组 成 ,因此 地 震 作 用 下 两 种 等 效 茜 性 
阻尼 模型 的 计算 结果 近似 相等 ,上 且 近 似 等 于 结构 阻 
尼 模 型 中 剔除 发 散 项 后 满足 物理 意义 的 计算 结 
进而 从 理论 上 说 明了 方法 的 合理 性 。 


图 区 
< ”一 一 一 损耗 因子 0.10 
pp 
1.10 A 一 -一 损耗 因子 0.20 
公信 公公 人 从 从 全 全 全 全 
， “损耗 因子 0.40 


动力 放大 系数 比值 


SO 1.0 10.0 
频率 比 
(a) VET 


一 一 一 损耗 因子 0.10 

1.10 四 于 020 
机 

，- - <- - 损耗 因子 0.40 


全 全 全 从 全 全 全 全 


202307.00057V1 


动力 放大 系数 比值 


01 1.0 10.0 
频率 比 
(b) VET 

图 2 不 同等 效 阻尼 模型 的 动力 放大 系数 比值 


.2 The variation laws of dynamic magnification factors 


ThinaXiv 


based on different equivalent damping models 

以 损耗 因子 为 0.10, 自 振 频率 为 10 rad/s 的 单 
自由 度 体 系 为 例 ,分 别 采用 VFT、VET、 基 于 结构 阻 
尼 模 型 的 时 域 计 算 方法 (complex time-domain damp- 
ing model ,CT) 和 基于 结构 阻尼 模型 的 频 域 计算 方法 
(complex frequency-domain damping model, CF ) 计算 
单 自由 度 体系 在 El Centro 波 、 天 津 波 作用 下 的 时 域 
动力 响应 时 程 ,计算 结果 见 图 3 和 图 4。 其 中 ,VFT 
的 阻尼 系数 为 0.998 8 ,VET 的 阻尼 系数 为 1. 000 0。 
如 图 3 和 图 4 所 示 , 随 着 地 震 作 用 时 间 的 增加 ,CT 
的 时 域 计算 结果 逐渐 出 现 发 散 现象 ,原因 如 第 1 节 
所 述 结构 阻尼 模型 的 自由 振动 通 解 中 存在 指数 增加 
项 ;VFT、VET 和 CF 的 时 域 动力 响应 时 程 一 直 保 持 
稳定 收敛 , 且 计 算 结 果 近 似 相等 (图 3 ~4) ,对 应 的 
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响应 峰值 均 小 于 5% ( 表 1) ,进而 证 明了 VFT 和 
VET 的 正确 性 。 以 上 计算 结果 表明 ,与 结构 阻尼 模 
型 的 时 域 计 算 方 法 相 比 , 基 于 等 效 和 共性 阻尼 模型 的 
计算 方法 可 有 效 解 决 时 域 发 散 问 题 。 


位 移 /cm 


时 间 /s 
(a) 位 移 时 程 


[> 
b= 
让 
[= 


一 CT ---VvFT 


Ri 


oY 1000 
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9 
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其 
网 
兵 -1 000 

-2 000 

2 4 6 8 10 12 14 16 
时 间 /s 
(b) 加 速度 时 程 


图 3 ”也 Centro 波 作用 下 的 结构 动力 响应 时 程 
Fig.3 Time-history of structural dynamic response 


under El Centro wave 


4 
人 
3 
2 
E 中 
半 0 
HE 
蕊 -1 
3 
3 
-4 
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400 
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南 
其 -200 
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时 间 /s 
(b) 加 速度 时 程 
到 4 ”天津 波 作用 下 的 结构 动力 响应 时 程 


Fig.4 Time-histories of structural dynamic 


responses under Tianjin wave 
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表 1 不 同方 法 计算 结构 响应 峰值 的 对 比 
Tab.1 Comparisons of peak structural responses with different methods 
El Centro 波 天 津 波 
类 别 

CF VFT VET CF VFT VET 
位 移 峰 值 /cm 4.7874 4.7285 4.7272 3.3066 3.3206 3.3200 

相对 误差 /% 三 1.23 1.26 = 0.42 0.41 
加 速度 峰值 / (cm . s-?) 583.1070 578.4930 578.3730 400.5690 398.3170 398.1880 

相对 误差 /多 = 0.79 0.81 一 0.56 0.59 


4.2 多 自由 度 体系 算 例 


针对 非 比 例 阻尼 体系 ,基于 阻尼 衰减 等 效 的 
Rayleigh 阻尼 模型 和 基于 阻尼 耗 能 等 效 的 Rayleigh 
阻尼 模型 的 计算 流程 图 如 图 5 所 示 , 依 据 两 种 等 效 
阻 兹 模型 可 计算 结构 的 时 域 动力 响应 。 


材料 损耗 


Rayleigh 阻 


基于 阻尼 衰减 等 效 的 
Rayleigh 阻 尼 模 型 


对 域 动力 响应 


图 5 等 效 Rayleigh 阻尼 模型 的 计算 流程 图 
Fig.5 The flow chart of equivalent Rayleigh damping mode 

以 6 自由 度 体系 为 例 ,质量 和 刚度 分 布 如 图 6 
所 示 , 模 型 A 中 上 部 2 层 结构 的 材料 损耗 因子 为 
0.04, 下 部 4 层 结构 的 材料 损耗 因子 为 0. 10 ;在 模型 
A 的 基础 上 ,保持 质量 刚度 不 变 , 仅 增 大 材料 的 阻尼 
性 能 ,可 进一步 得 到 模型 B, 上 部 2 层 结构 的 材料 损 
耗 因 子 为 0.10, 下 部 4 层 结构 的 材料 损耗 因子 为 
0.20。 然 后 ,分 别 采 用 基于 阻尼 衰减 等 效 的 Rayleigh 
阻尼 模型 时 域 计 算 方法 (Rayleigh time-domain damp- 
ing model based on free vibration response ,RFT) \ 基 于 
阻尼 耗 能 等 效 的 Rayleigh 阻尼 模型 时 域 计 算 方法 
(Rayleigh time-domain based on the equivalent damp- 
ing energy dissipation , RET) 和 CF 计算 模型 A 和 模 
型 B 在 了 Centro 波 、 天 津 波 作 用 下 结构 顶层 的 动力 


基于 阻尼 衰减 等 效 的 
Rayleigh 阻 尼 模 型 


— 上 时 域 动力 响应 
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响应 时 程 ,结果 如 图 7 ~ 10 所 示 。 在 计算 过 程 中 ,可 
直接 选取 结构 的 前 两 阶 振 型 自 振 频 率 确定 Rayleigh 
阻尼 和 矩阵 ( 算 例 为 简单 结构 ) ,RFT 和 RET 对 应 的 阻 
尼 系 数 如 表 2 所 示 。 


2 000 kg 


1.0X10 Nm 
1.5X10 Nm 
2.0X10 Nm 
2.5X10 Nm 
3.0X10 Nm 
”| 3.5X10 Nm 


图 6 多 自由 度 模型 示意 图 


Fig.6 Schematic diagrams of the multi-degree of freedom models 
表 2 不 同等 效 Rayleigh 阻尼 模型 的 阻尼 系数 


Tab.2 Damping coefficients of different equivalent 


Rayleigh damping models 


_ RFT RET 
类 别 

a B a B 
模型 A 0.1847 0.0063 0.1849 0.0063 
模型 B 0.3625 0.0145 0.3641 0.0146 


El Centro 波 作 用 下 ,RFT、RET 和 CF 对 应 的 结 
构 顶 层 动力 响应 时 程 近 似 相 等 (图 7 ~ 图 8)。 如 表 
3 所 示 ,对 于 模型 A,RFT 的 顶层 位 移 峰 值 加 速度 峰 
值 的 相对 误差 分 别 为 1. 02% .0.70% ;RET 的 顶层 
位 移 峰 值 、 加 速度 峰值 相对 误差 分 别 为 1.05% 、 
0.71% 。 对 于 模型 B, RFT 和 RET 的 顶层 位 移 (加 
速度 ) 峰值 相对 误差 分 别 为 1. 18% (0. 94% )、 
1. 30% (1.19% ) 。 

对 于 天 津 波 作用 下 的 模型 A, 如 图 9 所 示 , 在 
0 ~8 s 时 ,RFT、RET 和 CF 的 顶层 位 移 时 程 和 加 速度 
时 程 差 异 较 大 ,分 析 其 原因 ,RFT 和 RET 为 时 域 计 
算 方法 ,对 应 的 计算 结果 包含 瞬 态 响应 和 稳 态 响应 ， 
CF 为 频 域 计 算 方 法 ,对 应 的 计算 结果 仅 包含 稳 态 响 


] 力 学 学 报 


isyiv 信 人 作 划 : 
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应 。 在 8s 以 后 , 瞬 态 啊 应 逐渐 衰减 消失 ,RFT、RET ”顶层 位 移 ( 加 速度 ) 峰值 相对 误差 为 2. 66% 
和 CF 的 动力 响应 时 程 近似 相等 。 对 于 天 津 波 作用 ”0.09% )。 对 于 模型 B,RFT 的 顶层 位 移 ( 加 速度 ) 
下 的 模型 B,RFT、RET 和 CF 的 动力 啊 应 时 程 近似 。 ”峰值 相对 误差 为 0.78% (2. 14% ) ;RET 的 顶层 位 移 
相等 。 如 表 4 所 示 , 对 于 模型 A, RFT 的 顶层 位 移 。 (加速度) 峰值 相 对 误差 为 0.93% (2.35% ) 。 
( 加速度) 峰值 相 对 误差 为 2.65% (0.10% ) ;RET 的 

表 3 ”也 Centro 波 作用 下 不 同方 法 计算 的 结构 响应 峰值 


Tab.3 Peak structural responses with different methods under El Centro wave 


模型 A 模型 B 
类 别 
CF RFT RET CF RFT RET 
位 移 峰 值 /em 39.947 5 39.540 6 39.528 6 34.7867 34.375 5 34.335 0 
相对 误差 /% 村 1.02 1.05 一 1.18 1.30 
加 速度 峰值 / (cm . s-2) 478.453 6 475.1130 475.0400 370.0540 366.5740 365.6620 
相对 误差 /% 一 0.70 0.71 一 0.94 1. 19 


> 表 4 天津 波 作用 下 不 同方 法 计算 的 结构 响应 峰值 
i> Tab.4 Peak structural responses with different methods under Tianjin wave 
© 模型 A 模型 B 
© 类 别 
© CF RFT RET CF RFT RET 
~ 位 移 峰 值 /em 23.8408 23.2093 23.2069 20.0542 19.898 1 19.8678 
蝗 相对 误差 /% 一 2.65 2.66 一 0.78 0.93 
(af 
CW 速度 峰值 A( em . s 2?) 437.037 4 437.4700 437.4261 361.1605 353.4467 352.6651 
©O 相对 误差 /% 一 0. 10 0.09 一 2.14 2.35 
§ 
; 泥 
全 
- 
© 
0 10 20 30 40 
时 间 /s 时 间 /s 
和 (Ga) 位 移 时 程 
600 — RFT 
立 S 
训 让 
昌 站 
0 10 20 30 40 
时 间 /s 时 间 /s 
(b) 加 速度 时 程 (b) 加 速度 时 程 
图 7 El Centro 波 作用 下 模型 A 的 顶层 动力 响应 时 程 到 8 El Centro 波 作 用 下 模型 B 的 顶层 动力 响应 时 程 
Fig.7 Time-histories of top floor dynamic responses of Fig.8 Time-histories of top floor dynamic responses of 
Model A under El Centro wave Model B under El Centro wave 
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(b) 加 速度 时 程 
后 图 9 天 津 波 作 用 下 模型 A 的 顶层 动力 响应 时 程 


Sy 
CBig.9 Time-histories of top floor dynamic responses of 
ss 


Model A under Tianjin wave 
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人 时 间 /s 


(a) 位 移 时 程 


/(cm* s) 


加 速度 


时 间 /s 
(b) 加 速度 时 程 


图 10 天 津 波 作用 下 模型 B 的 顶层 动力 响应 时 程 


Fig. 10 Time-histories of top floor dynamic responses 


of Model B under Tianjin wave 

由 于 结构 阻尼 模型 存在 自由 振动 通 解 发 散 、 有 
阻尼 自 振 频率 随 损耗 因子 增 大 而 增 大 等 理论 缺陷 ， 
导致 基于 阻尼 衰减 的 计算 结果 和 基于 阻尼 耗 能 的 计 
算 结 果 存 在 一 定 差异 ,两 种 结果 不 能 形成 完全 统一 
的 形式 。 但 两 种 等 效 阻尼 模型 的 计算 结果 近似 相 
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等 ,与 自由 衰减 和 耗 散 能 量 均 能 反映 结构 阻尼 耗 能 
本 质 的 理论 是 一 致 的 。 

综 上 所 述 , 由 于 结构 阻尼 模型 的 频 域 计算 可 视 
为 准确 值 , 上 述 计算 结果 对 比 证 明了 等 效 Rayleigh 
阻尼 模型 时 域 计算 方法 的 正确 性 。 但 是 ,基于 结构 
阻尼 模型 的 频 域 计算 方法 仅 适用 于 初始 时 刻 处 于 静 
止 状 态 ( 初 值 条 件 为 零 ) 的 结构 ,等 效 Rayleigh 阻尼 
模型 作为 一 种 时 域 阻尼 模型 ,可 考虑 瞬 态 啊 应 的 影 
响 ,计算 初始 时 刻 处 于 非 项 止 状态 的 结构 动力 响应 ， 
进一步 拓展 了 结构 阻尼 模型 的 适用 范围 。 

作为 非 比 例 阻 尼 结 构 体 系 的 近似 时 域 计算 方 
法 ,等 效 Rayleigh 阻尼 模型 可 适用 于 非 比例 阻尼 特 
性 较 弱 的 混合 结构 。 与 基于 阻尼 衰减 等 效 的 Ray- 
leigh 阻尼 模型 相 比 ,基于 阻尼 耗 能 等 效 的 Rayleigh 
阻尼 模型 计算 结果 近似 相等 , 且 不 涉及 求解 复 特征 
值 ,计算 过 程 更 简单 ,建议 优先 选用 。 


5 结 论 


结合 黏 性 阻尼 模型 和 结构 阻尼 模型 的 优点 ,本 
研究 提出 了 与 结构 阻尼 模型 等 效 的 Rayleigh 阻尼 模 
型 ,并 通过 算 例 分 析 得 出 以 下 结论 。 

1) 分 别 依据 阻尼 衰减 和 阻尼 耗 能 ,构建 了 等 效 
于 结构 阻尼 模型 的 单 自由 度 体系 黏 性 阻尼 模型 。 分 
析 表 明 ,等 效 昔 性 阻尼 模型 可 克服 结构 阻尼 模型 时 
域 计 算 结果 发 散 、 初 值 条 件 不 易 确定 等 缺陷 。 

2) 在 单 自 由 度 体 系 等 效 黏 性 阻尼 模型 的 基础 
上 ,进一步 构建 了 多 自由 体系 的 等 效 Rayleigh 阻尼 ， 
保留 了 结构 阻尼 模型 中 阻尼 和 矩阵 易 构造 的 优点 ,可 
由 材料 损耗 因子 、 结 构 质量 矩 阵 和 刚度 矩阵 确定 对 
应 的 阻尼 系数 ,避免 了 模 态 阻尼 比 的 测定 ,计算 过 程 
简便 ,可 直接 用 于 计算 非 比 例 阻尼 体系 的 时 域 动力 
响应 。 

3 ) 结 构 阻 尼 模 型 仅 适用 于 频 域 分 析 , 等 效 Ray- 
leigh 阻尼 模型 的 时 域 计算 方法 可 考虑 瞬 态 响应 , 进 
一 步 拓 展 了 结构 阻尼 模型 的 适用 范围 。 与 基于 阻尼 
衰减 等 效 的 Rayleigh 阻尼 模型 相 比 ,基于 阻尼 耗 能 
等 效 的 Rayleigh 阻尼 模型 计算 结果 近似 相等 , 且 不 
涉及 复 模 态 分 析 , 计 算 过 程 更 简单 。 
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